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Abstract. We present a novel description of nuclear many-body systems, both for nuclear matter and finite nuclei, 
emphasizing the connection with the condensate structure of the QCD ground state and spontaneous chiral symmetry 
breaking. Lorentz scalar and vector mean-fields are introduced in accordance with QCD sum rules. Nuclear binding 
arises from pionic fluctuations, using in-medium chiral perturbation theory up to three-loop order. Ground state proper- 
ties of 16 and 40 Ca are calculated. The built-in QCD constraints reduce the number of input parameters significantly 
in comparison with purely phenomenological relativistic mean-field approaches. 
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^ • 1 Introduction 



The description of nuclear many -body dynamics must ultimately 
be constrained by the underlying theory of the strong interac- 
tion - Quantum Chromodynamics (QCD). Previous phenomeno- 
logical steps with this goal in mind have been taken by Quan- 
tum Hadrodynamics (QHD) 1 1 1. In the mean field (Hartree) ap- 
proximation, such an approach is equivalent to a model with 
local four-point interactions between nucleons 121314 1. Models 
based on QHD have been successfully applied to describe a va- 
riety of nuclear phenomena over the whole periodic table, from 
light nuclei to superheavy elements (see Ref. [5| for a recent 
review, and references therein). 

While this phenomenological success is impressive, an un- 
derstanding of its foundations in QCD is still missing. The 
multitude of input parameters in QHD models is usually not 
constrained by QCD considerations. Explicit pionic degrees of 
freedom are absent in most QHD type calculations, whereas 
it is obvious that pions, as Goldstone bosons of spontaneously 
broken chiral symmetry, must play an important role in the nu- 
clear many-body problem. Two-pion exchange effects are sup- 
posedly incorporated as part of the strong scalar-isoscalar field 
of QHD models, but in an ad-hoc manner without detailed ref- 
erence to the underlying irirNN dynamics. 

A general low-energy effective Lagrangian for nuclear sys- 
tems can been written down as a Taylor series in point cou- 
plings involving nucleon currents and their derivatives 16171 . A 
large number of coefficients must be determined in such an ef- 
fective field theory. The empirical data set of nuclear bulk and 
single-particle properties can be used to fix no more than six 
or seven of these parameters. Our approach is similar in spirit 
but proceeds with a different strategy, imposing as many QCD 



constraints as possible in order to minimize the number of free 
parameters. 

The success of relativistic mean-field phenomenology has 
been attributed primarily to large Lorentz scalar and vector nu- 
cleon self-energies |7|. There is evidence, in particular from 
nuclear matter saturation and from spin-orbit splittings in finite 
nuclei, that the magnitudes of these scalar and vector potentials 
are of the order of several hundred MeV in the nuclear interior. 
Investigations based on QCD sum rules |8 9 10 1 have shown 
how such large scalar and vector nucleon self-energies arise in 
finite-density QCD, at least qualitatively, through changes in 
the quark condensate and the quark density. Such QCD sum 
rule constraints will be one of the basic elements of our discus- 
sion. 

The second important ingredient is chiral pion-nucleon dy- 
namics. In Ref. 1 1 1 1 the equation of state of isospin-symmetric 
nuclear matter has been calculated using in-medium chiral per- 
turbation theory. At nuclear matter saturation density, the Fermi 
momentum kf and the pion mass represent comparable 
scales, and therefore pions must be included as explicit de- 
grees of freedom in the description of nuclear many-body dy- 
namics. The calculations have been performed to three-loop or- 
der and incorporate the one-pion exchange Fock term, iterated 
one-pion exchange and irreducible two-pion exchange. The re- 
sulting nuclear matter equation of state is expressed as an ex- 
pansion in powers of the Fermi momentum kf. The expansion 
coefficients are functions of kf/m^, the dimensionless ratio 
of the two relevant scales. The calculation involves one single 
parameter, the momentum space cutoff A which encodes NN- 
dynamics at short distances. With A ~ 0.65 GeV adjusted to 
the energy per particle E(kfo) = —15.3 MeV, the calculated 
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equation of state gives the density po — 0.178 fm -3 , the com- 
pression modulus K = 255 MeV, and the asymmetry energy 
A(kf ) = 33.8MeV at saturation. 

Based on these observations, our "minimal" approach for 
nuclear matter and finite nuclei starts from the following hy- 
pothesis: 

A) The nuclear matter ground state is characterized by large 
scalar and vector nucleon self-energies of approximately equal 
magnitude and opposite sign, arising from the in-medium change 
of the scalar quark condensate and the quark density. 

B) Nuclear binding and saturation result from chiral (pio- 
nic) fluctuations superimposed on the condensate background 
fields. These pionic fluctuations are calculated according to the 
rules of in-medium chiral perturbation theory. 

As concerns hypothesis A), finite-density QCD sum rules 
I8I9I1 0I predict the scalar and vector potentials to be each about 
300 — 400 MeV in magnitude at nuclear matter saturation den- 
sity po. The same QCD sum rule analysis, taken to leading or- 
der, also suggests the ratio of scalar to vector fields to be close 
to —1. We shall argue that this is indeed a valid starting point, 
though not yet capable of producing nuclear binding. Hypoth- 
esis B) asserts that binding and saturation is ruled primarily 
by explicit tttt exchange dynamics based on known properties 
of irN interactions and calculable using systematic methods 
of chiral effective field theory - at least as long as the Fermi 
momentum hp is small compared to the characteristic scale, 
4^/^ ~ 1 GeV, associated with spontaneous chiral symmetry 
breaking in QCD. 

Our aim in this paper is thus to study the interplay between 
condensate background fields and perturbative chiral fluctua- 
tions, both rooted in the spontaneous symmetry breaking pat- 
tern of QCD, in forming nuclei. We will demonstrate that this 
scenario works at large, once a single scale parameter is set to 
reproduce nuclear matter at equilibrium. Whereas in first ap- 
proximation the condensate potentials do not play a role for 
the saturation mechanism, we will show that they are essential 
for the description of ground states of finite nuclei. We restrict 
ourselves here to gross features of isospin-symmetric (N = Z) 
nuclei and relegate further fine tunings as well as the N > Z 
case to forthcoming work. 



2 Model for nuclear matter and finite nuclei 

2.1 Lagrangian 

Our approach is defined by the following (isoscalar) Lagrangian, 
relevant for N = Z nuclei: 

+2 G .s(p) V>V>V>V> - -G v {p)$"f^-4)^4> 

+^D s (p)d^d^-^D v (p)d^j^d v ^^ 



(1) 



Here, ip is the nucleon spinor field, M is the (free) nucleon 
mass and the subscripts S and V refer to the scalar and vector 
type interactions, respectively. The vector potential and field 



strength tensor of the electromagnetic field are denoted and 
F^ h '. The coupling parameters of the four-nucleon contact in- 
teractions and the derivative terms are assumed to be functions 
of the nucleon density p. These coupling strengths include con- 
tributions from condensate background fields and pionic (chi- 
ral) fluctuations, to be specified. We will formally work at the 
mean field level using Eq.(l), with the understanding that fluc- 
tuations beyond mean field are encoded in the density depen- 
dent coupling strengths. 

The single-nucleon Dirac equation derived from the La- 
grangian Eq. Q by variation with respect to ip, reads: 

[ 7 „ {id» -H*- E£) - (M + S s + So,)] 1> = , (2) 

with the nucleon self-energies defined by the following rela- 
tions: 



E 1 * = G v f - D v Uf - e A» 

Z s = -Gs(W>)+DsU(iH>) 
dD s , 



1+T 3 



Er 



dp 



(3) 
(4) 
(5) 



1 dGy .,, . 1 dDy . . . . 



dD 



v 



dp 



{d v j a )u a {d»f) . 



(6) 



where = Tpj^ip is the nucleon current, and the velocity 
it M is defined by pu^ = j M . In addition to the usual vector 
IJV and scalar S s self-energies, the density dependence of the 
vertex functions Gs(p), Gy(p), Ds(p) and Dy(p), produces 
the rearrangement contributions £r s and 1121 . The inclu- 
sion of the rearrangement self-energies is essential for energy- 
momentum conservation and the thermodynamical consistency 
of the model 1121131 . 

The ground state of a nucleus with A nucleons is the prod- 
uct of the lowest occupied single-nucleon self-consistent sta- 
tionary solutions of the Dirac equation Eq.(|2j. The ground state 
energy is the sum of the single-nucleon energies plus a func- 
tional of the scalar density 

A 

ps = ^2 i>ki>k 
fc=i 

and of the nucleon (vector) density 

A 

p = Yl ^ k 

k=l 

calculated in the no-sea approximation, i.e. the sum runs only 
over occupied positive-energy single-nucleon states with wave 
functions tpk . 



2.2 Nuclear matter 

The energy density £ and the pressure P of isospin symmet- 
ric nuclear matter are calculated from the energy-momentum 
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tensor as 



£=£o 



GspI 



-Gyp 2 



(7) 



P 



d£ 



£ = fi*p — £ + -Gyp 



; G sPs 



ldG s a 



ldGv 3 



(8) 



It should be pointed out that, while these expressions are 
formally derived in the mean field (Hartree) approximation from 
the Lagrangian (1), they incorporate exchange effects and fluc- 
tuations beyond mean field. In particular, the pionic fluctua- 
tions to be described in more detail in section 2.4 are calculated 
at three-loop order which includes Fock terms from one-pion 
exchange as well as all possible exchange terms related to two- 
pion exchange. These effects are transcribed into the density 
dependence of the couplings Gs.v{p)- 

The free quasi-particle contribution is given by 



£o 



(27T) 



d 3 ky/k 2 +M* 2 = h$n*p + M*p s ), 



\k\<kf 

with the effective chemical potential 



and the effective nucleon mass 

M* = M - G s ps ■ 



(9) 



(10) 



(ID 



The baryon density is related to the Fermi momentum kf in 
the usual way, p — 2k 3 /3ir 2 , and the expression for the scalar 
density reads 



PS 



(2tt) 3 



d A k- 



M* 



M* 



\k\<k f 

k f p* - M* 2 In 



Vfc 2 + M* 2 
kf+p* 



M* 



(12) 



Note that, in contrast to the energy density, rearrangement con- 
tributions appear explicitly in the expression for the pressure. 
The general form of the vertex functions Gs(p) and Gy{p) 

is 



Gs(p) = 4° } - AGs(p) 
G v {p) = Gf + AG v (p) 



(13) 
(14) 



where Ggy are terms governed by the QCD condensates, and 
AGs,v (p) refer to the pionic fluctuations, re-expressed as den- 
sity dependent corrections to the mean-fields. 

2.3 Constraints from QCD condensates 

The in-medium QCD sum rules relate the changes in the scalar 
quark condensate and the quark density due to the finite baryon 



density with the scalar and vector self energies of a nucleon in 
the nuclear medium. In leading order, which should be valid 
below and around nuclear matter saturation density, one finds 
for these condensate parts of the nucleon self energies 1 9 1 : 



.(0) 



— H< QQ >p 



A 



< 11 >^ 



87T 2 



B 



A 2 



E 



(0) 



64tt 2 



32n 2 
~A 2 ~ 



p 



PS 

m d 
(15) 



(16) 



where Ab ~ 1 GeV is a characteristic scale, the Borel mass, 
entering in the QCD sum rule analysis. For typical values of 
the nucleon sigma term cat and the current quark masses m„ 
and md, the ratio 



S 



(0) 



a N 



.(0) 



A{m n 



m d 



(17) 



is close to — 1 (take, for example, — 45 MeV and m u + 
md — 12 MeV), with uncertainties at the 20% level. Using the 
Hellmann-Feynman theorem in combination with PCAC (the 
Gell-Mann-Oakes-Renner relation) to derive the in-medium scalar 
quark condensate, one finds in the Fermi gas approximation 1 8 1 : 



,(0) 



M* -M 



9 £9 PS 



which implies: 



Ss\p ) ^ -350 MeV 



cat Ps 
50 MeV 



or, identifying Eg 



G s ps- 



G 



(0) 



11 fm 2 



<7 N 



50 MeV 



at ps — Po = 0.16 fm 



(18) 



(19) 



(20) 



Up to this point our discussion is based on the leading terms 
of the QCD sum rule, the ones involving the dimension-four 
condensates (m q qq) and (G^ V G^ V ). The density dependence 
of the gluon condensate is in fact weak and need not be con- 
sidered. The influence of higher-dimensional condensates has 
been discussed in great detail in Ref. |9|. Uncertainties arise 
primarily from contributions of four-quark condensates. It is 
common practice to approximate those four-quark condensates 
assuming factorization which introduces potentially large and 
uncontrolled errors. We can refrain from this discussion be- 
cause our explicit treatment of scalar irir fluctuations removes 
at least part of these uncertainties. 



2.4 Pionic (chiral) fluctuations 

This brings us next to the constraints from chiral pion-nucleon 
dynamics on the density dependent parts of Eqs.(13,14). If we 
follow the assumption, also made implicitly in Ref. [11|, that 



in nuclear matter E 



(0) 



,(0) 



v 5 — —^ v at p = po, the density de- 
pendent couplings of the pionic fluctuation terms AGs{p) and 
AGy (p) are determined by equating the corresponding self- 
energies in the single-nucleon Dirac equation (0 with those 
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calculated using in-medium chiral perturbation theory (CHPT) 
inRef.fll]: 



AG s ( P )ps = S^ FT (kf,p) , 



(21) 



I d AGs 2 
AG v (p)p + -—p s 



1 OAGy 2 _ v^CHPT 



(kf,p) 



(22) 

We have indicated here that the £™ PT (p, p) depend explic- 
itly on the nucleon momentum p. 

The energy per particle E(kf) in nuclear matter gives, via 
the Hugenholtz-van Hove theorem, the sum of the scalar and 



_ ^CHPT 



S 



(kf,p) 



vector nucleon self-energies U(kf,kf) 
Sy FT (kf, p) at the Fermi surface p — kf up to two-loop 
order, as generated by chiral one- and two-pion exchange 1 14 1. 
The difference Sg HPT (p, p) — Ey PT (p, p) is calculated from 
the same pion-exchange diagrams via the anti-nucleon single 
particle potential in nuclear matter. Following a procedure sim- 
ilar to the determination of the nucleon-meson vertices of rel- 
ativistic mean-field models from Dirac-Brueckner calculations 
1151 . we neglect the momentum dependence of Sg } v PT (p, p) 



and take their values at the Fermi surface p — kf. A polyno 
mial fit up to order k^ is performed, and the self-energies art 
then reexpressed in terms of baryon density p = 2/c 3 /37r 2 : 



S 



CHPT 



(k f ,p) = (c 50 + csi p 1/3 + c S2 p z/3 ) P (23) 



££ HPT (fc/, p) = (cvo + cvi p l/i + c V 2 P 11 ' 6 ) P - (24) 

The values of the coefficients are: c$o = —2.805 fm 2 , csi — 
2.738 fm 3 , c S 2 = 1-346 fm 4 , c VQ = -2.718 fm 2 , c vi = 
2.841 fm 3 , and Cy2 = 1-325 fm 4 . The resulting expressions 
for the density dependent couplings of the pionic fluctuation 
terms are 



AG s (p) = c 50 + csi p 1/3 + c S 2 p 



o 2 /3 



(25) 



AG V (p) = Cy 



-(6cvi 



■c S i)p 1/3 + \(3c V2 -c S 2) P 2/3 ■ 

(26) 

In deriving the expressions for AGs(p) and AGy (p) we have 
set ps ~ p on the left hand sides of of Eqs. Oil and \22\ . 
Although the relation between the scalar and baryon density 
depends on the Fermi momentum, this approximation is jus- 
tified for p < po- With the density dependent couplings d25t 
and i26\ . Eqs. to H2\ produce a nuclear matter equation of 
state which is very close to the one calculated successfully in 
CHPT. This is shown in Fig. ;1; where we compare the nuclear 
matter equation of state calculated from Eqs. and (jSJi, with 
the one obtained using in-medium CHPT II II . Correspond- 
ing ground state properties: the binding energy per particle, the 
saturation density, the compressibility modulus, and the asym- 
metry energy at saturation, are compared in Tabled Small dif- 
ferences arise mainly because the momentum dependence of 
the CHPT self-energies has been frozen in Eqs. (23, 24). This is 
a well known problem which has been extensively discussed, 
for instance, in Ref. 1 15 1. By fine tuning just two of the pa- 
rameters in Eqs. i25\ and (126b we could, of course, reproduce 
the CHPT equation of state of Ref. II II exactly. In the present 
work, however, we prefer not to perform any such tuning of 
parameters. 



2.5 Finite nuclei 

Having constrained Gg y from QCD condensates, and hav- 
ing adjusted one short-distance scale parameter appearing in 
the pionic fluctuation couplings AGs(p) and AGv{p) to the 
equation of state of isospin symmetric nuclear matter, we pro- 
ceed to calculate finite nuclei. In this work we only consider 
the isoscalar channel and calculate the ground states of le O 
and 40 Ca. In addition to Gs(p) and Gy(p), two new quantities 
appear specifically for finite nuclei, namely the couplings of 
the terms involving derivatives in the nucleon fields in Eq.Q: 
Ds{p) and Dy{p)- Guided by dimensional considerations, we 
introduce the ansatz 



Ds(p) = ^- and Dy(p) 



Gv(fi) 
A 2 



(27) 



where A is again a characteristic mass scale delineating short 
and long distance phenomena. In the present calculation we 
simply choose A w 0.65 GeV, the same value that has been 
used for the momentum space cutoff in the in-medium CHPT 
calculation of the nuclear matter equation of state. In this way, 
and we emphasize this point, no new parameters are needed in 
the calculation of finite nuclei. 

In the first step we have calculated the ground states of le O 
and 40 Ca with G s {p) = AG s {p) and G v {p) = AG v {p), i.e. 
we have set the couplings to the condensate background fields 
to zero. The nuclear dynamics is then completely determined 
by chiral (pionic) fluctuations. The interesting result is that the 
calculated total binding energies are within 5 — 8% of the ex- 
perimental values, but the resulting radii of the two nuclei are 
too small (by about 0.2 fm). This is because the spin-orbit part- 
ners (lp 3 /2, IP1/2) and (ld 5 / 2 , ld 3 / 2 ) are practically degener- 
ate: chiral two-pion exchange dynamics alone, although it can 
provide the attraction necessary to bind nuclei, does not pro- 
duce the proper spin-orbit interaction. This is shown in Fig. [2] 
where we display the neutron and proton single-particle levels 
in le O and 40 Ca, calculated in the limit G ( ^ v = 0. 

The spin-orbit degeneracy is removed by including the self- 
energies which arise from the changes in the scalar quark con- 
densate and quark density. In the second step we have adjusted 
the couplings gP and GP to the binding energies and the 
charge radii of l6 (D and 40 Ca. We emphasize again that, up to 
this point, our calculation of both the nuclear matter equation 
of state and the binding energies of finite nuclei includes only 
one adjustable parameter: A — 0.647 GeV. Even though Gg 

and GP were varied independently, the minimization proce- 
dure tends to favour cancellation of the contributions from the 
corresponding large scalar and vector self-energies. This hap- 
pens because there is already enough binding from pionic fluc- 
tuations, and therefore sP = — Ey ' represents a very good 
approximation for the condensate potentials. The final values 



Gp = 10.52 fm z and G^> = 10.00 fm" 2 should be compared 



,(0) 



with the estimate Eg. (1201 . and with the leading order coeffi- 
cients of the pionic (CHPT) terms cso and cyo- This is a re- 
markable result which indeed supports the "minimal scenario" 
with condensate background fields plus pionic fluctuations as a 
very reasonable starting point. 
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In Table|2]we compare the calculated binding energies and 
charge radii with the corresponding empirical values. The ab- 
solute deviations between theory and experiment are 7.8% and 
5.3% for the binding energies, and 2.5% and 4.6% for the charge 
radii of le O and 40 Ca, respectively. In Fig. [5] the neutron and 
proton single-particle levels in le O, calculated with the inclu- 
sion of the condensate potentials, are compared with the ex- 
perimental levels. We notice that this calculation reproduces 
about 2/3 of the empirical spin-orbit splitting. While this re- 
sult is clearly not yet satisfactory at a quantitative level, it also 
indicates the necessary steps for further improvements. The ap- 
proach is so far not fully self-consistent, in the following sense. 
By construction, the condensate potentials do not contribute to 
the binding of nuclear matter which is accounted for almost 
entirely by chiral two-pion exchange dynamics. This leaves no 
room for increasing the scalar (S) and vector (V) condensate 
background contributions such that, by their difference S — V, 
the large spin-orbit splitting can be reproduced quantitatively. 
The obvious solution is to treat the chiral (two-pion exchange) 
fluctuations and the condensate self-energies on the same level, 
both for nuclear matter and for finite nuclei. Furthermore, it is 
necessary to go beyond the leading order in the QCD sum rules 
and examine higher order density dependence for the conden- 
sate self-energies. These points will be considered in a forth- 
coming analysis of a generalized point-coupling model con- 
strained by QCD sum rules and in-medium CHPT 

We note in passing that our model for finite nuclei could, in 
principle, be based on an alternative chiral approach to isospin- 
symmetric nuclear matter proposed by Lutz et al [16|. In ad- 
dition to pion-exchange, their approach includes contributions 
from a zero-range NN-contact interaction treated beyond the 
mean-field approximation (i.e. the contact interaction is also it- 
erated with lir exchange). However, while the nuclear matter 
equation of state of Ref. 1161 is comparable to the one used 
in the present work, the single-nucleon potential resulting from 
that approach has several unrealistic features 11171 . For instance, 
the single-nucleon potential at zero momentum {7(0, fc/o) ~ 
—23 MeV is not sufficiently attractive, and the total single- 
nucleon energy does not increase monotonically with momen- 
tum (implying a negative effective mass). The corresponding 
scalar and vector self-energies resulting from this particular 
scheme would therefore not be suitable for applications to fi- 
nite nuclei. 



dimensional QCD condensates. The calculations are also ex- 
panded to cover a wider range of finite nuclei with extensions 
towards N > Z systems. 

This work has been supported in part by BMBF, GSI, DFG and INFN. 
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3 Conclusion 

The effective Lagrangian Eq. Q, with couplings governed by 
scales of low-energy QCD, gives a good description of both 
symmetric nuclear matter and finite N = Z nuclei at a level 
better than 10%, even without detailed fine tuning of param- 
eters. While nuclear binding and saturation are almost com- 
pletely generated by chiral (two-pion exchange) fluctuations, 
strong scalar and vector fields of equal magnitude and opposite 
sign, induced by changes of the QCD vacuum in the presence 
of baryonic matter, drive the large spin-orbit splitting in finite 
nuclei. Considering the constraints from QCD condensates and 
chiral dynamics that keep the number of adjustable parameters 
at minimum, our results are quite encouraging. Investigations 
are now being generalized to include corrections from higher 
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Table 1. Nuclear matter saturation properties: binding energy per nu- 
cleon, saturation density, incompressibility, and the asymmetry energy 
at saturation. The first row corresponds to the in-medium CHPT cal- 
culation including one- and two-pion exchange 1 11 1. The EOS dis- 
played in the second row is obtained when the resulting CHPT single 
nucleon potentials are mapped on the self-energies of the relativistic 
point-coupling model with density dependent couplings. 
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MODEL 


E/A (MeV) 


po (/m- 5 ) 


K (MeV) 


A (MeV) 


CHPT 1111 


-15.26 


0.178 


255 


33.8 


PC-DD 


-14.51 


0.175 


235 


36.6 



Table 2. Binding energies per nucleon E/A [MeV] and root mean 
square charge radii r c [fm] of ie O and 40 Ca. The experimental val- 
ues, shown in the second and third column, are compared with the 
results of the present calculation. 



E exp /A 


exp 


E/A 


r c 


lb O 7.98 


2.74 


8.60 


2.80 


40 Ca 8.55 


3.48 


8.10 


3.64 
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Fig. 2. Neutron and proton single-particle levels in 16 O and 40 Ca 
calculated in the relativistic point-coupling model. The calculation is 
performed using only the contribution from chiral one- and two-pion 
exchange to the density dependence of the coupling parameters (i.e. 
Gf v = 0inEqs.(13,14)). 
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Fig. 1. Binding energy per nucleon for symmetric nuclear matter as 
a function of the baryon density. The solid curve (KFW-2002) is the 
EOS calculated in Ref. 1111 by using in-medium CHPT. The EOS 
displayed by the dotted curve (PC-dd) is obtained when the resulting 
CHPT nucleon potentials are mapped on the self-energies of the rela- 
tivistic point-coupling model with density dependent couplings. 
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Fig. 3. The neutron and proton single-particle levels in O calculated 
in the relativistic point-coupling model, are shown in comparison with 
experimental levels. The calculation is performed by including both 
the contributions of chiral pion-nucleon exchange, and of the isoscalar 
condensate self-energies, to the density dependence of the coupling 
parameters. 



